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Abstract 

In this paper we study the phase diagram of two Ising planes coupled by a stan- 
dard spin-spin interaction with bond randomness in each plane. The whole phase 
diagram is analyzed by help of Monte Carlo simulations and field theory arguments. 



1 Introduction 

In order to understand the role of impurities and inhomogeneities in real physical systems, 
many statistical models with quenched randomness have been proposed. The effect of 
randomness on continuous phase transitions has thus been of great interest for many 
years. The main question is to determine whether the randomness leaves the critical 
properties of the pure system unchanged. One prediction concerning models with random 
bonds comes from the Harris criterion [0, which states that bond randomness changes 
the values of the critical exponents only if the specific heat exponent a is positive. The 
effects of bond randomness have been studied intensively in the 2D Ising model and more 
recently in the g < 4 states Potts model. In the former case, the log singularity of the 
specific heat is transformed in a log{log) singularity 0, whereas in the latter case a new 
critical point with new critical exponents is reached 0. Furthermore, it has been recently 
realized that a weak bond randomness can also have strong effects on 2D spin systems 
possessing a first order transition. Following earlier work of Imry and Wortis [Q, Hui 
and Berker have shown, using a phenomenological renormalization group, that bond 
randomness (even infinitesimal) induces a second order transition (i.e. vanishing latent 
heat) in a system that would have undergone a first order one. This result has been proved 
rigorously by Aizenman and Wehr At this level, a question naturally arises concerning 
the universality class of this second order transition. In order to test these predictions, 
the q states Potts model (g > 4 is necessary to get a first order transition) with bond 
randomness has been recently extensively analyzed numerically 0], the conclusion being 
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that the transition is of second order type due to disorder and the critical behaviour 
depends on q. Analytical results on the q states Potts model are, as usual, much more 
difficult to obtain. However, this is not the case for those systems presenting what 
Cardy called a "weak fluctuation-driven ffist order transition" [^. The typical example 
Cardy provided is a system of critical Ising models coupled by their energy density. 
Without disorder, this perturbation drives the system in a strong coupling regime where 
perturbative renormalization group (RG) analysis fails. Yet, for this model, it can be 
shown non-perturbatively that this runaway flow is indeed associated to a massive regime 
(with a finite correlation length depending on the coupling). This is a way to define this 
"weak fluctuation-driven first order transition" phenomenon. The addition of weak bond 
randomness, even infinitesimal, makes the system flow to decoupled critical Ising 
models [§]. This result has been extended to the case of A^ coupled Potts models [P, |T0[] . 
The result was a non-Ising like second order transition, the universality class depending on 
the sign of the coupling between the models. This class of examples is interesting because 
some infinitesimal bond randomness is able to change the "weak fluctuation- driven ffist 
order transition" into a continuous one. Concretely speaking, the RG runaway flow is 
completely turned over by disorder preventing a strong coupling regime. Nevertheless a 
lot of questions still remain open. What happens if we take a different coupling than 
the energy density one? Does infinitesimal disorder will also change the order of the 
transition? 

To give some insights to these questions, we would like to study the following system: 
two Ising planes coupled by a spin-spin interaction where disorder (dilution or bond- 
randomness) is added in each plane. When the two Ising planes are critical, the spin- 
spin perturbation which is much stronger that the energy-energy perturbation drives the 
system in a strong coupling regime. It has been shown using integrable Toda field theories 



that a finite correlation length is generated |]TT| , and the full spectrum of the theory can be 
computed too. In this sense, the situation is analogous to the model treated by Cardy ||], 
except that now the coupling is different and much stronger. Therefore, this model can 
bring some clues to the aforementioned questions. Moreover, it is worth noticing that A^ 
Ising planes (with A^ ^ 1) coupled by the spin-spin interaction is a way to reach a 3D 
Ising model. N = 2 may be regarded as a kind of ffist step in this direction. Therefore, 
we expect Aizenman-Wehr theorem not to apply here. 

The plan of the paper is the following. In section 2, we present analytical field theory 
arguments showing that a finite amount of disorder is necessary to change the "weakly 
fluctuation-driven order ffist order transition" into a continuous one. In section 3, we 
conffim it by Monte Carlo numerical simulations. In section 4, we will summarize the 
results and discuss the whole phase diagram of the model. 
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2 Analytical results 



2.1 Pure case 

In this section, we first summarize results available in the pure case. The Hamiltonian 
then reads 



3 ' « J 



where a]'"^ = ±1, the first sum is over the nearest neighbours on a 2D square lattice, 
the same coupling J has been chosen for the two models and a coupling J' has been 
considered between the two planes. This model can be described in the continuum limit 
by the following action 

A = A], + Aj. + m J d'^x{ei{x) + e2ix)) + g J d^x (Ji{x)a2ix) , (2) 

where A]^ denotes the action of the critical Ising model, ei{x) the Ising energy operator, 
(Tj the Ising spin operator, m (x {T — T^) and g oc J'. 

When g = 0, the action correspond to two decoupled massive Ising models which can 
be described by free massive Majorana fermions or by free Dirac fermions or by a sine- 
Gordon model. If the strongly relevant coupling g is switched on at the m = critical 
point, the system is known to be driven into a massive regime whose mass spectrum has 
been computed exactly |]TT|. When g ^ 0, the bosonic representation is more suitable (in 
terms of Majorana fermions representation, this interaction is indeed non-local), and the 
generalized sine-Gordon reads 

A = J d^x^iVcpf + mcos{p^) + gcos{a(l)±5) , (3) 

with l3 = v^ivr, a = y/n and (5 = f ■ 

In order to obtain this expression, we write the Hamiltonian in terms of Majorana 



fermions and then use standard bosonization formula (see for example |12] for a review). 



It is worth noticing that the addition of a four spin interaction like K J2 



<i,3> 

the original Hamiltonian (|I]) (which therefore defines a generalized Ashkin- Teller model) 
would not change the form of the action (H) except the substitution (v^ivr, ^tF) — ^ 
(/?, \P) with /? = /?(ir). 

The action (|^) is a priori hard to study since it contains two strongly relevant per- 
turbations and it is unfortunately non integrable. Notice that g is much stronger that m 
and scales around the critical point like g ~ |m| ''''^ ~ |T — Tc\''^^. A perturbed conformal 
field theory approach is unappropriate and beyond feasibility because the flow is quickly 
driven in a strong coupling regime. 
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Figure 1: Evolution of the potential when varying g at fixed mass m. Criticality is 
reached for a critical value Qc- 



The action (0) which is one particular case of the double frequency sine-Gordon model 
has been qualitatively studied in the general case by Delfino and Mussardo |13|. It has 
been argued by looking at the solitonic structure in both weak and strong coupling regime 
(m 3> (7 and m ^ g) that when the ratio ^ = ^ is a rational number and 5 = ^, a phase 
transition should occur for m > {i.e. T > Tc). Concerning our particular case, we can 
easily convince ourselves about this transition by looking at the evolution of the shape of 
the potential V = — mcos(/?0) — gsm{f3(j)/2) which passes from a periodic two degenerate 
minima to a periodic one absolute minimum situation when increasing the ratio g/m (see 
Figure [^). From a Ginzburg-Landau point of view, it suggests that the transition is in 
the universality class of namely is of Ising type. 

Let us come back to the original Hamiltonian (|l|). In the limit g oo, the spins aj 
and af are locked to the same value and the model reduces to only one Ising model with 
a critical temperature 2Jc. For a generic value of g the same scenario should occur for 
a critical temperature Jc < Tc{g) < 2Jc since g tends to order the system in contrast to 
temperature. Therefore, there should be a critical line in the (T, g) plane interpolating 
between the fixed points (J^O) and (2Jc, oo). In the g ^ J limit, we can rewrite the 
Hamiltonian (|l|) in the basis at = aj; Ti = (jjaf and integrate over the r variable using 
a mean field approximation (we refer to |T^ for a more detailed analysis in a different 
context). It remains one Ising model with approximate effective coupling 2J(1 — c^) ~ 
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2Jefj, with c a positive number satisfying a self- consistent equation. The Ising criticahty 
is reached when Jeg = Jc- This analysis is valid only close to the (2 J^, oo) Ising fixed 
point, namely in the strong coupling limit. 

According to Zamolodchikov c-theorem |[T^ (without disorder the theory is still uni- 
tary), the flow goes from the ( Jc? 0) fixed point with central charge c = 1 toward the 
ending Ising fixed point (2 Jc, cxd) with central charge c = |. From the continuum limit 
point of view, it is worth noticing that the elementary excitations around this fixed point 
have nothing to do with the ones around the c = 1 fixed point (the correspondence is 
indeed non local explaining why a standard perturbation theory around the c = 1 fixed 
point is hopeless). 



2.2 Disordered case 

Let us now add bond randomness or dilution in each Ising planes. In the two limits (7 = 
and g = oo the problem reduces to the well known disordered Ising model and there are 
no new fixed points. Can we say something in the general situation? The phase diagram 
is parametrized by three variables, T, g and A the strength of the disorder. In the 
planes A = and g = 0, the shape of the transition separating the ferromagnetic phase 
from the paramagnetic one is known respectively from results concerning the disordered 
Ising model and from the pure case analysis described above. Therefore, by a continuity 
arguments, we expect a critical surface joining the two Ising critical lines. Let us focus on 
the T = Tc plane. In the continuum limit, the effect of disorder is to change m into m{x) 
in the action (^ with m(a;) = 0, 'm{x)m{x') = A6{x — x') for a Gaussian disorder. This 
action is very hard to study with usual tools and we can only give qualitative arguments. 
At fixed value for g, when m{x) > m*{g) = g/A, it implies a local high temperature 
region surrounded by low temperature regions. In the limit A ^ m* (g) (still at fixed g) , 
the phenomenology of the critical disordered Ising model is recovered with a mixture of 
low and high temperature regions. This suggests that there should exist a finite critical 
value A*{g) of the disorder strength able to disorder the system. In fact, the situation 
resembles the one of a disordered Ising model with m{x) = tuq > 0. This is clearly a 
non-perturbative phenomenon. We conjecture this transition to be of Ising type governed 
by the Ising fixed point (with c = |) at (2 Jc, oo, A = 0) with a logilog) behavior for the 
specific heat due to the disorder. Another non trivial fixed point at A* > cannot 
a priori be excluded but is very unprobable according to the symmetries of the model 



under consideration. Note that contrary to the examples treated in p, ^, [T^ where some 
infinitesimal disorder changes the weakly driven first order transition into a continuous 
one, here a finite value of the disorder is needed! 
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Figure 2: Susceptibility data rescaled with pure 2D Ising model exponents: z/ = 1 and 
7 = 7/4. In the left inset we report the values for the critical coupling Qc versus the 
disorder amount A, together with the best power law fit. In the right inset we show, as a 
function of the g values where the specific heat (filled symbols) and the susceptibility 
(empty symbols) have a peak. Here A = 0.25 and the Binder cumulants crossing point 
is gc = 0.743 (filled square). 

3 Numerical analysis 

To put on firmer grounds these qualitative arguments, we have performed Monte Carlo 
simulations at T = Tc. Let us present our numerical results. We have considered two 2D 
Ising models, with bond dilution A in both planes, i.e. P{J) = A6{J) + (1 — A)5( J — 1). 
The disorder is independent in each plane. We let the planes interact via a spin- spin 
coupling term, whose strength is g. 

We have performed numerical simulations with 3 different lattice sizes {L = 20, 40, 60) 
and for three values of the disorder: A = 0.05,0.15,0.25. Lattice sizes may seem small 
with respect to previous studies on the disordered Ising model [|16[, however they seem 
to be in the scaling regime (see below) and so the extrapolation to large sizes is safe, 
at least up to the precision we are interested in. For any disorder value we have found 
clear evidences that, increasing the value of the coupling g, the system undergoes a phase 
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transition from a paramagnetic to a ferromagnetic phase. The critical couphng values 
can be well fitted by a power law: gd^) oc A^'^ (see left inset in Fig. 0, where the errors 
are smaller than the symbol size). Note that the value of this exponent is very close to 
the theoretical one obtained by scaling arguments close to the c = 1 fixed point, which 
is 7/4. 

For every A, the critical coupling has been determined as the crossing point of 
the Binder cumulants. This crossing point has almost no finite size effects. Two more 
estimations of gc can be obtained from the g values where the specific heat and the 
susceptibility have a peak. In fact these values converge to gc in the large L limit. However 
they show much larger finite size effects than the Binder cumulants crossing point. Then, 
in order to have a better estimation of gc, we have used the data for every lattice size and 
we have considered finite size corrections through the formula gdL) = gc — AL~^^'^, with 
V = 1. We have found that our data are compatible with that formula for any disorder 
value and the gc estimation coincides with the one extracted from the Binder cumulants 
(see right inset in Fig. |^). The effect of increasing A simply reflects in a larger value for 
the non-universal coefficient A. 

In order to verify that the transition is governed by the pure 2D Ising model fixed 
point, we have analyzed our numerical data using finite size scaling and the critical 
exponents of the pure 2D Ising model, i.e. u = 1 and 7 = 7/4. We have found very good 
data collapse for any disorder value A. In particular in Fig. |^ we show the collapse of 
the rescaled susceptibility for the larger disorder value A = 0.25. For smaller disorder 
amplitudes the scaling is even better and the scaling function is more peaked signaling 
that the critical region is narrower. 

4 Discussion and conclusions 

In this paper we have presented analytical and numerical evidences that in two coupled 2D 
Ising models the "weak fluctuation-driven flrst order transition" changes to a continuous 
one only when a flnite amount of disorder is added. Moreover, we have seen that even for 
flnite and large disorder amplitudes, the phase transition is in the 2D disordered Ising 
universality class (Ising transition modifled by logarithm corrections). The results we get 
can be summarized in the schematic 3D phase diagram depicted in Fig. |^. Our analysis 
strongly suggests that almost all the critical surface shown in Fig. ^ is governed by the 
2D Ising flxed point, except the curve in the g = plane (governed by the c = 1 flxed 
point) and the one in the T = plane (governed by a percolation flxed pointQ). Because 

Note that the T = critical hne is discontinuous at g = 0, as can been seen in Fig. H In fact for 
every g ^ the bond percolation threshold is l/V^, while for g = is 1/2 
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Figure 3: Phase diagram in the variables T,g,A, where A is the bond dilution. The 
critical surface separates the ferromagnetic region (below) from the paramagnetic one 
(above). The arrows indicate the flow direction. We have numerically studied the system 
along the line T = shown in the figure. 

of the known flow lines on this surface (the ones with an arrow in Fig. we believe 
that any transition across the surface is governed by the c = 1/2 fixed point located 
in (2 Jc, oo) (the rightmost big dot in Fig. ||) because the c = 1 fixed point is strongly 



repulsive in the g direction (even if it is attractive on the g = plane [0). For an even 
clearer numerical check, a numerical estimation of the central charge would be wellcome. 

As a conclusion, we would like to mention that several directions remain open. First 
of all, a theoretical framework to describe quantitatively such non-pertubative transition 
is needed. Secondly it would be nice extending this analysis for > 2 coupled Ising 
models and finally it is worth performing a similar analysis for correlated disorder. 
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